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ABSTRACT 

We present phase-resolved spectroscopy and extensive survey photometry of the detached 
eclipsing binary system XYCet, which is composed of two metallic-lined stars. We measure 
their masses to be 1.773 ±0.016 and 1.615 ±0.014 M Q and their radii to be 1.873 ±0.035 and 
1.773±0.029 R , resulting in logarithmic surface gravities of 4.142±0.016 and 4.149±0.014 
(cgs). We determine effective temperatures of 7870 ± 115 and 7620 ± 125 K. The projected 
rotational velocities are 34.4 ± 0.4 and 34.1 ± 0.4kms _1 , which are close to synchronous. 
Theoretical models cannot match all of these properties, but come closest for a solar helium 
and metal abundance and an age in the region of 850 Myr. We obtain the individual spectra 
of the two stars by the spectral disentangling method, and compare them to synthetic spectra 
calculated for the measured effective temperatures and a solar chemical composition. Both 
stars show enhanced abundances of iron-group elements and clear deficiencies of Ca I and 
Sc II, confirming their classification as Am stars. We also find strong overabundances of Zrll 
and the rare-earth species La II, Cell and Ndll, a hallmark of chemically peculiar A stars. 
XY Cet is a prime candidate for detailed spectroscopic analyses of metallic-lined stars whose 
masses and radii are known to accuracies of 1-2%. 

Key words: stars: binaries: eclipsing - stars: fundamental parameters - stars: binaries: spec- 
troscopic - stars: early-type — stars: chemically peculiar 



1 INTRODUCTION 

Eclipsing binary star systems are of immense value as our pri- 
mary source of empirical measureme nts of the properties of stars 
dAndersenl Il99ll : iTorres et al.l |2010). The masses and radii of 
those systems which show double-lined spectra can be measured 
using purely geo metrical arguments and to precisions of bet- 
ter than 1% (e.g. IClausen et al.l I2009T) . With effective tempera- 
ture determinations they become excellent distan ce indicators both 
within the Milky Wa y and to external galaxies (Rib as etHi1l2005l ; 
Bonan oTet alj|200r3) . Detached eclipsing binaries (dEBs) are the 
pick of the bunch as their properties can be used to challenge stel- 
lar theory and single-star evolutionary models. 

One use of dEBs is to probe chemical peculiarities, such as the 
Am phenomenon in stars ( hitus & Morgan|[l94ol : IContil 1 1970h by 
measuring their physical properties to high accuracy. Am stars are 
A stars which show photospheric abundance anomalies t hought to 
be caused by radiative diffusion and gravitat ional settling dMichaud 
ll97d : lTurcotte et al.l200ollTalon et al.l2006h . These effects are able 
to operate in the radiative atmospheres possessed by A stars with 
rotational velocities slower than roug hly 100 km s" 1 dAbt & Levvl 



1 19851 : iBudail Il996l. fl997» . Am stars are preferentially found in 
short-period binaries where tidal effects have been abl e to slow 
their rotation dAbtll 1 96 ll [l965l: iCarquillat & Prieurll2007l) . 

This in turn means that Am stars are str ongly represented 
among the well-studied dEBsjj, such as /?Aur ( Southworth et al. 
|2007|) . V3 64Lac dTorres et alT[l999T) , V6 24Her dPoroerl 1 19841) . 
V459 Ca s dLacv et al.l |2004|), W W Cam dLacv et alj |2002|) and 
RRLyn dTomkin & Fekell l200d) . The prevailing viewpoint from 
these studies is that the Am phenomenon is a surface disease - the 
properties of Am stars can in general be matched by theoretical 
predictions just as well as those of normal A stars - although a high 
bulk metal abundance ( 5% by mass) was found fo r the component 
stars of WW Aurigae bv lSouthworth et al.ld2005ch . 

XY Ceti is a dEB containing Am stars which has been studied 
on many occasions but whose physical properties were not firmly 
established. In this work we present new spectroscopy and exten- 
sive photometric measurements from the SuperWASP survey from 
which we measure the physical properties of the component stars to 
high precision. This study highlights the many possibilities opening 
up for the study of eclipsing binaries of all types from the increasing 
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1 A catalogue of well-studied detached eclipsing binaries is maintained at 
http : // www .astro. keele. ac. uk/~ jkt/debdata/debs . html 
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Table 1. Identifica tions, location, a n d phot ometric indices for XYCet. 
References: (1) IPerrvman et al l Jl997l); ( 2) lH0g et alj Jl997l) ; (3) 
Cannon & Pickering! Jl918l) : (4) lArgelanderl dl903h : (5) This work; (6) 



Skrutskie et 



alJ toOrj) . 





XY Ceti 


Reference 


Hipparcos number 


HIP 13937 


1 


Tycho designation 


TYC 51-832-1 


2 


Henry Draper number 


HD 18597 


3 


Bonner Durchmusterung 


BD +02°460 


4 


RA (J2000) 


02 59 33.5333 


1 


Dec. (J2000) 


+03 31 03.274 


1 


Hipparcos parallax (mas) 


1.72 ± 1.52 


1 


Spectral type 


A8 Vm + F0 Vm 


5 




9.094 ± 0.024 


2 


V T 


8.800 ± 0.026 


2 


«^2MASS 


8.212 ± 0.023 


6 


^2MASS 


8.194 ±0.031 


6 


^2MASS 


8.115 ±0.021 


6 



number of long-term variability surveys, in particular those which 
have the photometric precision to search for extrasolar planetary 
transits. 



1.1 XY Ceti 

The objec t BD +02°460 was discovere d to be an eclipsing binary 
system bv lStrohmeier & Kniggej dl96ll) . who found that an orbital 
period of 1.390356 d satisfied their photographic observations but 
yielded a light curve with no dete ctable secondary minimum. The 
first photoelectric study of XY Cet dMorrison & Morrisonll968h re- 
vised the orbital period to 2.78071 2d, resulting in a light curve 
containing rather similar primary and secondary minima. Morri- 
son & Morrison analys ed their light curve using the method of 
iRussell&Merrillld 19521) . finding that the two stars were very sim- 
ilar. 

IPoppeJ d 19691) announced that both components of XYCet 
had a metallic- lined nature. This was independently discovered by 
Dr. A. Young (Popper Hl97ll) . The first - an d so far only - time- 
resolved spectroscopic study was presented by |Poppei1dl97ll) based 
on photographic material. Popper combi ned his spectroscopic 
results with the photometric elements of iMorrisori & Morrison 
dl968l) to make the first determination of the masses and radii of the 
component stars of XY Cet, finding the two objects to have similar 
properties (masses 1.75 and 1.63 Mq and radii 1.95 and 1.89 R©). 
Popper gave a rough spectral classification of A2m + FOm. 

However, an analys i s of n ew photometric data for XYCet 
led ISrivastava & Padalial dl975) to radii w hich were quite dif- 
ferent (2.14 and 1.61 Rm ubkazakil dl978l) revisited the photo- 
metric data obtained by b oth [Morrison & Morrisonl d 19681) and 
ISrivastava & Padalial JT975) and resolved this discrepancy in favour 
of Popper's results (radii 1.96 and 1.78 Rq), in the process finding 
that XY Cet has a small but significant orbital eccentricity. How- 
ever, a substantial dispa rity between the rad ii of the two compo- 
nents was resurrected by Ramella et al. ( 1980), through an analysis 
of the data obtained bv lSrivastava & Padalial dl975h . These authors 
used, for the first time in the case of XY Cet, 'modern light curve 
modelling techniques' to find radii of 2.10 and 1.74 Rq. These 
values we re supported by ISrivastaval dl987l) . who reanalysed the 
data from ISrivastava & Padalial d 19751) to find values of 2.13 and 
1.76 R . 



Since the apparent resolution of varied measurements of the 
radii of the component stars of XY Cet in favour of a notable dif- 
ference between the t wo objects, there h as been only a little interest 
in this binary system. ISrivastaval J 1988b presented a period study 
of XY Cet, finding five changes in period but providing no clear 
evidence for the ir existence. Finally, o ne spectrum of the object 
was obtained bv lGlazunova et af. (2008), who measured the rota- 
tional velocities of both components. Their observation was taken 
at an orbital phase of 0.957, where the spectral lines of the two 
components would be partially overlapping, so it is not surprising 
that their results are inconsistent with the values we find in Sect.[6] 
The basic observational properties of XY Cet are given in TableQ] 
Throughout this work we refer to the primary star as star A and the 
secondary star as starB. Star A is hotter, larger and more massive 
than starB. 



2 DATA ACQUISITION 
2.1 Spectroscopy 

Spectroscopic observations were obtained in 2002 October using 
the 2.5 m Isaac Newton Telescope (INT) on La Palma. The 500 mm 
camera of the Intermediate Dispersion Spectrograph (IDS) was 
equipped with a holographic 2400 lines mm -1 grating. An EEV 
4k x 2k CCD was used and exposure times were 300 s. From mea- 
surements of the full width at half maximum (FWHM) of arc lines 
taken for wavelength calibration we estimate that the resolution 
is 0.2 A (1.8 pixels). A total of 38 spectra were taken covering 
the interval 4230-4500 A, with estimated signal to noise ratios of 
roughly 150 per pixel. One spectrum was taken covering 4710- 
4970 A, aimed at the H/? line. 

The r eduction of a ll spectra was undertaken using optimal 
extraction dHornel 1 19861) as implem ented in the software tools 



PAMELA and MOLLY 



So) as impleme 
-lMarsl] dl989l) . 



2.2 Photometry 

Our main photometric datas et comes from the W ASP survey for 
transiting extrasolar planets dPollacco et al.l l2006) which operates 
two instruments, one in La Palma and one in South Africa. Each 
instrument consists of eight small telescopes on a common mount, 
and each telescope comprises a Canon 200 mm tele photo lens, 
a custom filtejf] resembling g+r dFukugita et"aIH l996), and an 
E2V 2kx2k CCD. Both instruments studied the area surrounding 
XY Cet extensively in the 2008 and 2009 o bserving seasons. The 
data were reduced by a dedicated pipeline dPollacco et al.ll2006T) 
which performs aperture photometry in three softwa re apertures 
and t hen applies the SysRem detrending algorithm (Tam uz et al.l 
120051) to the photometry from the second of these apertures. 

The WASP observations of XYCet consist of 15 227 data- 
points in total, of which 13 702 were observed in good conditions 
during the 2008 and 2009 seasons. From preliminary light curve 
fits (see Sect.[5} we established that the four different photomet- 
ric aperture/detrending approaches gave results in very good agree- 
ment with each other and that the SYSREM-detrended data had the 
lowest scatter. We therefore adopted the detrended data for all fur- 
ther analyses. These data were iteratively clipped to remove outliers 

2 PAMELA and MOLLY were written by Dr. Tom Marsh and are available 

at |htt p : / /www . Warw ick . ac . uk/go/trmarsh| 

3 Seewww.superwasp.org 
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Table 2. Literature times of minimum light of XY Cet and the observed 
minus calculated (O — C) values of the data compared to the ephemeris 

derived in this work. 

Refere nces: (1) Strohmeier & Knig gel ll96ll); (2) iMorrison & Morrisonl 
19681) : Ol lBaldwirj ll973l): W lSrivastaval ll987l): (5) Kurpinska-Winiarska 
199611: (6) lLocheJ 1 19723): (7) lLocheil il972bl): (8) iDiethelrri jl974a| 



(9) lDiethelnJ ll974bl) ; (l oTlAl-Naimiv etalj ll978l): q i )|Paschd ll989; 
?T ll989bl) : (13)lDie 



(12)|PascheM "l989b): ( 1 3) lDiethelnJ 11 1 99 lb : (14) lBlattleil 1 19921) : (15l lZeidal 
12004F (16) This work. 



Time of minimum 


Cycle 


D — L 


Reference 




number 


(riJIJ ) 




26734.271 ±0.04 


-4185.5 


0.00724 


1 


26945.596 ± 0.04 


-4109.5 


-0.00206 


1 


26945.611 ±0.04 


-4109.5 


0.01294 


1 


27365.467 ± 0.04 


-3958.5 


-0.01895 


1 


27386.362 ± 0.04 


-3951.0 


0.02069 


1 


27397.423 ± 0.04 


-3947.0 


-0.04117 


1 


28127.413 ±0.04 


-3684.5 


0.01127 


1 


36843.514 ±0.04 


-550.0 


-0.03740 


1 


36943.559 ± 0.04 


-514.0 


-0.09812 


1 


36943.604 ± 0.04 


-514.0 


-0.05312 


1 


36850.510 ±0.04 


-547.5 


0.00681 


1 


36850.556 ± 0.04 


-547.5 


0.05281 


1 


36903.379 ± 0.04 


-528.5 


0.04224 


1 


37316.292 ±0.04 


-380.0 


0.01913 


1 


38372.949 ± 0.005 


0.0 


0.00461 


2 


39878.701 ± 0.01 


541.5 


-0.00030 


3 


40529.387 ± 0.005 


775.5 


-0.00150 


4 


40532.161 ±0.005 


776.5 


-0.00821 


4 


40543.298 ± 0.005 


780.5 


0.00593 


4 


40557.194 ±0.005 


785.5 


-0.00164 


4 


40889.495 ± 0.01 


905.0 


0.00397 


5 


40906.177 ±0.005 


911.0 


0.00168 


4 


40917.301 ± 0.005 


915.0 


0.00283 


4 


40931.204 ±0.005 


920.0 


0.00225 


4 


40967.357 ± 0.01 


933.0 


0.00597 


5 


40981.270 ±0.01 


938.0 


0.01539 


5 


41348.300 ± 0.01 


1070.0 


-0.00892 


6 


41587.445 ± 0.01 


1156.0 


-0.00537 


7 


42071.284 ±0.01 


1330.0 


-0.01070 


8 


42289.565 ± 0.01 


1408.5 


-0.01579 


9 


43453.3049 ± 0.005 


1827.0 


-0.00492 


10 


44166.546 ± 0.01 


2083.5 


-0.01709 


5 


46762.359 ± 0.01 


3017.0 


-0.00110 


5 


46773.496 ± 0.01 


3021.0 


0.01305 


5 


47471.435 ±0.01 


3272.0 


-0.00730 


11 


47535.380 ± 0.01 


3295.0 


-0.01873 


12 


48233.375 ± 0.005 


3546.0 


0.01692 


13 


48561.483 ±0.005 


3664.0 


0.00061 


5 


48625.4378 ± 0.005 


3687.0 


-0.00103 


14 


52279.2978 ± 0.0032 


5001.0 


0.00009 


15 


52949.4487 ± 0.0025 


5242.0 


-0.00121 


15 


54787.49993 ± 0.00015 


5903.0 


0.00014 


16 


55157.33430 ± 0.00026 


6036.0 


-0.00044 


16 



at the 4<7 level, resulting in a final tally of 9111 good datapoints 
from the 2008 season and 4091 from the 2009 season. 

To provide individual UBV colours and an independent check 
on the WASP data we acquired the photoelectric U BV light curves 
of XY Cet obtained bv lSrivastava & Padalial Jl975l) . The individual 
data are unfortunately not tabulated in that work, but are plotted as a 
function of orbital phase (their Fig. 1). We have manually extracted 
the data from this plot by displaying it on a computer screen and 



mouse-clicking on each datapoint. This procedure was not straight- 
forward because the model light curve is plotted as filled circles 
on the same axes as the individual observations are plotted as open 
circles. This has the twin effects of improving the overall appear- 
ance of the light curves without having to increase the informa- 
tion content of the data, and of obscuring some observations with 
the symbols representing the model light curve. It should also be 
borne in mind that any imperfections in defining the axis positions 
will cause the datapoints to be systematically shifted or stretched 
in phase or magnitude. A final caveat is that the values of the dat- 
apoints will be affected by discretisation due to the finite number 
of pixels on the computer screen, but this effect should be small 
(the plotting symbols were roughly five pixels wide). This proce- 
dure resulted in light curves containing 270 U, 263 B and 280 V 
datapoints; these numbers are rather smaller than the original work 
(335, 343 and 343, respectively) due to the difficulty of discerning 
overlapping plot symbols. The original data were made available to 
lOkazakil ll978h . who published normal points which can be used 
as an alternative to the digitised version obtained above. 



3 PERIOD DETERMINATION 

All available times of minimum light of XYCet were collected 
from the literature, and the orbital ephemeris given in the General 
Catalogue of Variable Star^ was used to determine the prelimi- 
nary cycle number of each minimum. Very few of the minimum 
timings are accompanied by an errorbar, so uncertainties were as- 
signed based on the size of the residuals for each type of observing 
equipment used. Two additional times of minimum were calculated 
from the combined WASP data, separated into the 2008 and 2009 
seasons, in the course of the analysis presented in Sect.[5] 

A straight line was fitted to the resulting cycle numbers 
and times of minima (TableO by % 2 minimisation. The resulting 
ephemeris is: 

Mini = HJD 2438372.9455(16) + 2.78071392(28) x E 

where the bracketed quantity is the uncertainty in the final digit 
of the preceding number. All uncertainties quoted in this work are 
standard errors unless otherwise stated. The residuals of the fit 
are plotted in Fig.[T]and give no indication of any form of period 
change. The above ephemeris is in good agree ment with that give n 
in the Eclipsing Binaries Minima Databas^ ( Kreiner et al ] l200lh . 
and shows no substantive evidence for departures from a constant 
period. There is also no evidence for the secondary minima being 
shifted from phase 0.5, implying that eccentricity is negligible. 



4 RADIAL VELOCITY ANALYSIS 

Our time-resolved spectroscopic dataset consists of 38 observations 
covering the wavelength range 4230-4500 A, of which six were 
taken during secondary eclipse when the two stars have very similar 
velocities. We have measured radial velocities (RVs) from the re- 
maining 32 spectra, concentrating on the 4360-4490 A wavelength 
range which contains a multitude of spectral lines but avoids the 
very broad H7 feature. We have considered three different meth- 
ods of measuring the velocity amplitudes of the component stars 



http : / /www . sai . msu . su/ groups /cluster /gcvs/gcvs / 
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Figure 1. Residuals of the times of minimum light around the ephemeris which best fits the observed times of minima. The filled circles refer to primary 
minima and the open circles to secondary minima. 



Table 3. Velocity amplitudes for the components of XYCet found from four different analysis methods. The probable errors given bv lPoppej fl97lh have 
been multiplied by 1 .48 to transform them into standard errors. 



Analysis type 


Template spectrum 


Uncorrected 


Corrected 




Light ratio 






K A (kms- 1 ) 


K B (kms" 1 ) 


^(kms" 1 ) 


K B 


(kms" 1 ) 




ONECOR 


HD 24740 


108.23 ±0.27 


118.65 ±0.45 


108.24 ±0.23 


118 


65 ± 0.28 




ONECOR 


HD 905 


108.12 ±0.26 


118.66 ±0.48 


108.10±0.21 


118 


65 ±0.21 




ONECOR 


HD 32115 


108.10 ±0.28 


118.72 ±0.43 


108.15 ±0.22 


118 


72 ± 0.22 




ONECOR 


HD 37594 


108.02 ±0.26 


118.74 ±0.40 


107.92 ±0.21 


118 


74 ± 0.20 




ONECOR 


HD 39945 


107.86 ±0.26 


118.66 ±0.44 


107.86 ±0.23 


118 


65 ±0.21 




TODCOR 


HD 24740 


108.42 ± 0.22 


118.89 ±0.34 


108.42 ±0.22 


118 


89 ± 0.22 


0.72 ±0.11 


TODCOR 


HD905 


108.32 ±0.21 


118.95 ±0.34 


108.13 ±0.09 


118 


89 ± 0.20 


0.69 ±0.13 


TODCOR 


HD 32115 


108.33 ±0.23 


119.12±0.33 


108.31 ±0.22 


119. 14 ±0.26 


0.73 ±0.10 


TODCOR 


HD 37594 


108.12 ±0.22 


119.24 ±0.30 


108.12 ±0.22 


119.22 ±0.24 


0.72 ±0.10 


TODCOR 


HD 39945 


108.02 ±0.23 


118.94 ±0.33 


108.03 ±0.21 


118 


94 ± 0.25 


0.70 ±0.11 


FDBINARY 




108.06 ±0.30 


119.22 ±0.35 










Pot>Der(1971) 




109.6 ±0.6 


117.6± 1.0 











of XY Cet, and this redundancy allows consistency checks and the 
assignment of robust measurement uncertainties. 

A large number of standard stars were observed using 
the same observational setup as for our target star. Inspection 
of these yielded five which have a similar appearance to the 
spectra of the components of XY Cet: HD 39945 (spectral type 
A5 V),HD 32115 (A8 IV), HD 37594 (A8Vs),HD905 (FOIV)and 
HD 24740 (F2IV). These will be used as template spectra in the 
analyses below. 



4.1 onecor: standard cross-correlation 

Numerical cross-correlation dSimkinlll974l: iTonrv & Davislll979l) 
is a standard approach for measuring RVs from the spectra of ce- 
lestial obje cts. We used our own implem entation of this method 
(ONECOR; [Southwort h & Clausen] 120071) after binning all spec- 
tra onto a common logarithmic wavelength scale. The cross- 
correlation functions (CCFs) were interactively assigned weights 
based on two factors: signal to noise and the velocity separation of 
the components of XY Cet. These weights were fixed for all sub- 
sequent analyses, and we have verified that their precise values do 
not have a significant effect on the resulting RVs. 

Each template spectrum was cross-correlated against the spec- 
tra of XYCet, and the positions of the two peaks were mea- 
sured using quadratic interpolation. The resulting RVs were fitted 



with spectroscopic orbits using the SBO qj code. Orbits were fit- 
ted for the two star s separately (see So uthworth et al.l l2004al and 
Popper & Hillll99ll) . Fits including orbital eccentricity yielded val- 
ues which were small and not significantly different from zero, as 
expected from the period study (Sect.|3}, so our final results were 
calculated with eccentricity fixed at zero. The outcome of this anal- 
ysis was velocity amplitude measurements for the two components 
of XYCet (Ka and Kb) and for each of the five template spectra 
(Tabled}. 

RVs found from cross-correlation analyses are known to show 
slight biases due to effects such as line blending jPetrie & Andrews] 
Il966l) and individual spectral lines being Dop pler-shifted into or 
out of the consider ed wavelength range (e.g. iTorres et al.lll997l : 
IClausen et al.l l2008h . These biases can be measured and thus re- 
moved by constructing synthetic composite spectra with known 
RVs and measuring them in the same way as the observed spec- 
tra. We have performed this analysis using our five observed tem- 
plate s pectra and the method discussed by Southwort h & Clausenl 
(120071) . We find that removing these biases from the measured RVs 
results in velocity amplitudes whose values are almost unchanged 
but whose uncertainties are noticeably lower; the results are given 
in Tabled 



Spectroscopic Binary Orbit Program, written by P. B. Etzel, 

|http : / /mintaka . sdsu . edu/ faculty/ etzel/| 
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Figure 2. RVs and best-fitting spectroscopic orbit measured using TODCOR 
and HD 905 as a template star. No velocity corrections have been made. 
Filled circles represent RVs for star A and open circles those for star B. 



4.2 todcor: two-dimensional cross-correlation 



The TODCOR algorithm, introduced by IZucker & Mazehl ( 1 19941) 
calculates a two-dimensional CCF for a double-lined spectrum us- 
ing two template spectra. This approach is aimed at avoiding line 
blending when the two target stars have q uite different spectral 
characteristics, although South worth & Clau sen (2007) found that 
it is no better than ONECOR for spectroscopic binaries containing 
two similar stars. This is the case with XY Cet, but the relatively 
low rotational velocities of its components mean that line blending 
is not a significant problem. 

We used our own implementation of TODCOR 
dSouthworth et al.l l2004al) and the same method and template 
stars as for ONECOR in deriving Ka and Kb- This process 
included the measurement and removal of RV biases, using 
synthetic spectra constructed from the observed template spectra. 
In each case we used the same template spectrum for both target 
stars, as the two components of XYCet have similar spectral 
characteristics. 

One advantage of TODCOR over ONECOR is that the optimal 
light ratio between the two target stars can be calculated as part of 
the algorithm. These light ratios are given in Table|3]but come with 
a health warning. Both components of XY Cet show the Am phe- 
nomenon, which means their spectral line strengths are peculiar and 
do not reflect their overall chemical composition. The light ratios 
from TODCOR strictly refer to the overall line strengths rather than 
the integrated light of each star. We have also not adjusted them for 
the difference in equivalent widths caused by the different effective 
temperatures of the two stars. They therefore cannot be used as ob- 
servational constraints, but should be good enough to act as checks 
on the results of the light curve analysis below. 

The TODCOR results are in general slightly better than the 
ONECOR results: the Ka and Kb values are very similar but the 
bias correction has a smaller effect (Table|3}- For our final results 
we adopt the values calculated using HD905 as a template star. 
The corresponding RVs are given in Table|4] The uncorrected spec- 
toscopic orbits are plotted in Fig.|2] the RV corrections are shown 
in Fig.[3]and the final spectroscopic orbits are displayaed in Fig.|4] 



Table 4. Corrected RVs from a TODCOR analysis of the spectra of XYCet 
using as templates a spectrum of HD 905. The columns marked O — C give 
the observed minus calculated velocity residuals. 



HJD — 
2400 000 



Weight 



Radial velocities 
star A star B 



O-C 
A B 



52569.58310 


2 


-47.97 


90.50 


-0 


95 


-2 


90 


52569.64324 


1 


-36.55 


78.81 


-1 


67 


-1 


27 


52569.69624 


1 


-25.67 


66.71 


-2 


34 


-0 


70 


52570.46545 


2 


124.21 


-96.25 


-0 


63 


-1 


13 


52570.46911 


2 


127.24 


-94.01 


2 


19 


1 


34 


52570.47278 


2 


126.82 


-94.20 


1 


56 


1 


38 


52570.51954 


2 


127.92 


-96.77 





69 





97 


52570.52321 


2 


128.00 


-97.25 





66 





61 


52570.52686 


2 


128.08 


-97.72 





65 





24 


52570.55581 


2 


126.73 


-99.75 


-1 


21 


-1 


23 


52570.55948 


2 


127.34 


-99.15 


-0 


63 


-0 


60 


52570.56315 


2 


127.91 


-98.32 


-0 


09 





26 


52570.61455 


2 


126.81 


-98.74 


-0 


74 


-0 


65 


52570.61821 


2 


127.46 


-97.79 





00 





20 


52570.62188 


2 


126.99 


-98.25 


-0 


38 


-0 


36 


52570.66646 


2 


125.25 


-96.31 


-0 


37 


-0 


34 


52570.67013 


2 


124.63 


-96.06 


-0 


80 


-0 


30 


52570.67380 


2 


124.42 


-95.81 


-0 


81 


-0 


27 


52571.58690 


2 


-51.35 


97.69 





17 


-0 


65 


52571.59056 


2 


—52.21 


99.01 


_o 


02 


-0 


07 


52571.59423 


2 


-52.31 


99.11 





54 


-0 


70 


52571.65538 


2 


-63.19 


111.02 





01 


-0 


14 


52571.65904 


2 


-63.94 


111.77 


-0 


17 


-0 


01 


52571.66270 


2 


-64.29 


111.95 





05 


-0 


45 


52571.71929 


2 


-72.04 


122.68 





28 


1 


52 


52571.72295 


2 


-72.89 


122.14 


-0 


11 





47 


52571.72663 


2 


-73.54 


122.01 


-0 


29 


-0 


17 


52571.73029 


2 


-73.57 


123.49 





13 





82 


52571.73395 


2 


-73.79 


123.54 





35 





38 


52572.45708 


1 


-25.92 


73.45 


1 


83 


1 


20 


52572.46075 


1 


-25.91 


73.50 


1 


03 


2 


13 


52572.46442 


1 


-25.10 


72.65 


1 


03 


2 


17 




0.4 0.1 

Orbital phase 

Figure 3. RV corrections measured for the TODCOR analysis with HD 905 
as a template star. Filled circles represent the corrections for star A and open 
circles those for star B. 
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Figure 4. RVs and best-fitting spectroscopic orbit measured using TODCOR 
and HD 905 as a template star, including velocity corrections. Filled circles 
represent RVs for star A and open circles those for star B. 



Table 5. Final spectroscopic orbital parameters for XY Cet. 



Parameter Primary Secondary 

Velocity amplitude K( km s" 1 ) 108.3±0.4 118.9 ± 0.5 

Systemic velocity (km s _1 ) -6.0 ±1.3 -5.7 ±1.3 
Mass ratio q 0.9108 ± 0.0051 

asini(R ) 12.482 ± 0.035 

Msin 3 i(M Q ) 1.768 ±0.016 1.611 ±0.014 



4.3 Final spectroscopic results 

The full set of Ka and Kb measurements is collected in Table[5] 
and includes values measured by spectral disentangling in Sect.|6] 
The values demonstrate a good agreement between different anal- 
ysis methods and different template spectra. We al so find an ac- 
ceptable agreement with the results o f Poppeil (ll97ih . For our final 
values we adopt the Ka and Kb found via TODCOR with HD 905 
as a template, but we increase the uncertainties in these quanti- 
ties to account for the variations in velocity amplitudes throughout 
TablefJ] The systemic velocities are measured relative to HD905, 
for which we adopt a systemic velocity of —25.8 ± 1.3 km s - 
dNordstrom et al.| [2004>. Table[5]gives final orbital parameters. 



Table 6. Model parameters of the WASP light curve of XY Cet. 



Parameter 




Value 


Sum of the relative radii 




0.2918 ±0.0011 


Ratio of the radii 


k 


0.947 ± 0.032 


Central surface brightness ratio 


J 


0.806 ±0.043 


Orbital inclination (°) 


i 


87.66 ±0.12 


Linear LD coefficient star A 


u A 


0.49 ±0.10 


Linear LD coefficient star B 


f B 


0.40 ±0.12 


ecosoj 




0.00021 ±0.00022 


Relative radius of star A 


VA 


0.1499 ±0.0028 


Relative radius of star B 


rB 


0.1419 ±0.0023 


Light ratio 


(■b/Ia 


0.752 ±0.061 



Table 7. Combined and individual UBV and uvby indices for 
XYCet. The combined UB V and uvby measurem ents are taken from 
Poppe r & DumonJ i 19771) and lHilditch & Hill Il975h . respectively. 



Parameter 


Combined 


Star A 


StarB 


V 


8.78 ±0.07 


9.38 ±0.07 


9.70 ± 0.07 


B-V 


0.26 ±0.01 


0.21 ±0.02 


0.33 ±0.03 


U -B 


0.14 ±0.01 


0.16 ±0.02 


0.11 ±0.03 


V 


8.74 ±0.05 


9.27 ± 0.06 


9.77 ± 0.07 


b-y 


0.16 ±0.01 


0.14 ±0.04 


0.18 ±0.06 


nil 


0.23 ±0.01 


0.27 ± 0.05 


0.17 ±0.08 


Cl 


0.82 ±0.02 


0.87 ± 0.04 


0.76 ± 0.06 


Toff (K) 


7750 ± 240 


7960 ± 560 


7570 ± 730 


Spectral type 




A8 


F0 




-0.2 0.0 0.2 0.4 0.6 

Orbital phase 



Figure 6. Plot o f the UBV light curves of XYCet from 
ISrivastava & Padalial )l975h compared to JKTEBOP models in which 
only the surface brightness ratio is an adjusted parameter. 

5 LIGHT CURVE ANALYSIS 

We have modelled the WASP light curve of XY Cet using the JKTE- 

P7i j — i ir~t 

BOP codqjcode (Southworth et al. 2004a. b), which is based on the 
EBOP program dPopper & Etzelll98lt lEtzel 198li: lNelson & Davisl 
Il972h . The fractional radii of the two stars (ta = — and tb = 
— where a is the orbital semimajor axis) are reparameterised as 
their sum and ratio: ta + f b and k = — because these are more 
directly related to the morphology of the light curve. Aside from 
ta + tb and k we included the orbital inclination (i), and the ratio 
of the central surface brightnesses of the stars (J) as fitted parame- 
ters. The gravity darkening exponents were set to 1.0 and the mass 
ratio (which has a negligible effect on the photometric solution) 
was fixed to 0.91. Attempts to fit for third light turned up a small 
and sub-zero value, so this parameter was fixed at zero for the final 
solution. 



JKTEBOP is written in FORTRAN77 and the source code is available at 

|http : / /www .astro . keele . ac . uk/^ jkt / codes . html| 
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Figure 5. Plot of the WASP light curve of XY Cet (filled grey circles) compared to the best-fit model (black line). The residuals of the fit are plotted at the 
base of the figure with an offset from zero. 



For limb d arkening (LD) we used the linear law (e.g. 
ISouthwortit 2008) and the good quality of the light curve meant we 
were able to include the LD coefficient for both stars as fitted pa- 
rameters. An alternative solution using the quadratic law, with the 
linear coefficient for star included as a fitted parameter, produced 
results which were practically identical to those obtained via the 
linear law. The resulting linear LD coefficients, «a = 0.49 ± 0.10 
and ub = 0.40 ± 0.12, are a little lower t han the theore tically ex- 
pected value of 0.58 for the g+r passband dClaretl2004ab . This can- 
not be attributed to the metallic nature of the stars, as coefficients 
f° r [w] ~ +0.5 are very similar to those for solar metallicity. 

The orbital period was fixed to the value found in Sect.[3]but 
the time of primary mid-eclipse was allowed to vary. A check for 
orbital eccentricity returned a positive result, so the quantity e cos w 
(where e is eccentricity and uj is periastron longitude) was also in- 
cluded as an adjusted parameter. The final value for this quantity, 
ecosuj = 0.00021 ± 0.00020, is not significant and is too diminu- 
tive to affect either the photometric or spectroscopic results. The 
final solution of the WASP light curve is given in Tableland plot- 
ted in Fig.[5] 

To find the la uncertainties in the measure d parameters, we 
have r un both 10000 Monte Carlo simulatio ns dSouthworth et"al] 
l2004ch and a residual-permutation analysis (Jen kins et al.l 120021) . 
The latter approach is able to account for any correlated errors (red 
noise) in the WASP data. We find, as expected, that the residual- 
permutation analysis returns larger error estimates than the Monte 
Carlo simulations. The residual-permutation error estimates are 
therefore adopted as the final errorbars in Table|6] 



5.1 Checks, comparisons and colour indices 

We checked the JKTEBOP solution of th e WASP data using the 200 3 
version of the Wilson-Devinney code dWilson & Devinnevlll97lb . 
The data were phase-binned into 140 points with the eclipse phases 



sampled more finely than the other phases, the n fitted using the 
JKTW D wrapper to the Wilson-Devinney code tSouthworth et al.l 
1201 lb . The resulting parameter values are consistent with JKTEBOP 
to within the errors, indicating that both solutions are reliable. 

The ISrivastava & Padalial d 1975b light curves were modelled 
using JKTEBOP with the geometrical parameters (ta, »"b an d *) 
fixed, the LD coe f ficients set t o theoretically expected v alues 
l l Van Hammd 1 1 9931 : Iciaret 2000; Iciaret & Hauschildtll2003h . and 
only J and the out-of-eclipse brightness adjusted. A good fit was 
found, corroborating the WASP light curve solution, and was used 
to obtain the individual UBV magnitudes and colours of the two 
stars (Table|7J. These should be treated with caution due to the 
modest quality of the data and the way in which they we re ob- 
tained for this study. A plot of the lSrivastava & Padaliald 1975b light 
curves and the JKTEBOP constraine d fit is shown in Fig.[6] Analysis 
of the normal points published by Okazaki d 1978b yielded values 
which were very similar, ind icating that our digitised version of the 
ISrivastava & Padalial ([1975) data has some reliability. 

Stromgren uvby photometry dStromgrenl Il966b of XYCet 
was o btained at five different orbital phases by iHilditch & HilH 
dl975b . One of these is during primary eclipse, which allows in- 
dividual Stromgren indices to be obtained for the two stars. These 
Stromgren index measurements were converted into uvby light 
curves and fitted using JKTEBOP, with all parameters except J 
and the out-of-eclipse brightness fixed. From the resulting light 
ratios and the combined Stromgren indices we have obtained the 
indices for the individua l stars (Table|7J. Using the uvby grids of 
ISmallev&Kupk3dl997h we find effective temperatures for the two 
stars of T eBA = 7960 K and T effB = 7570 K From the calibration 
of b — y versus spectral type dStromgrenl 1 9631) we find types of A8 
and F0. 

As a final investigation of the light curve we calculated a pe- 
riodogram of the residuals of the best JKTEBOP fit using the Pe- 
RIOD04 package dLenz & B reger 2004i). There is a peak of height 
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Figure 7. Plot of the disentangled spectra of star A (blue points, offset by 
+0.5) and starB (green points) around the H7 line. The best-fitting syn- 
thetic spectra for the two stars are shown with thin red lines . 



5 mmag at 1 cycle per day, and lower peaks at the multiples of this 
frequency, which are certainly of instrumental origin. There are no 
other clear periodicities in the frequency range of 7 Doradus pul- 
sations or S Scuti pulsations, to limiting amplitudes of 2 mmag and 
1 mmag respectively. The components of XY Cet are therefore not 
found to be pulsating stars. 



6 ATMOSPHERIC PARAMETERS FROM SPECTRAL 
DISENTANGLING 

Our spectroscopic dataset is well suited to a spectral disentangling 
analysis, in which the individual spectra and velocity amplitudes 
of the two stars are extracted from a set of phase-reso lved spectra 
dSimon & Sterol 1994t|pavlovski & Hensbergell2009l) . The result- 
ing disentangled spectra contain the total signal of the input spec- 
tra and are well suited to derivation of th e atmospheric parameters 
and chemical abundances of the stars (e.g. lPavlovski & Southworth] 
l2009t|Pavlovski et all2009h . 

We carefully normalised and velocity-binned the s pectra and 
disen tangled them using the FDBINAR"i0 package Jlliiic et all 
|2004|) . This analysis was initially performed both with and with- 
out the six spectra taken during secondary eclipse. We found that 
the former approach resulted in much more tightly constrained so- 
lutions, as these spectra contain information on how much light is 
lost from the secondary star at two orbital phases during secondary 
eclipse. The velocity amplitudes we find are included in Table[3] 

We measured the projected rotational velocities of the two 
stars by fitting the di s entang led spectra and the UCLSYN program 
dSmallev et all 1200 it ISmifhl Il992h , finding v A sin i = 34.4 ± 
0.4 km s" 1 and UBsini = 34.1 ± 0.5 km s -1 . These are close to 
the equatorial rotational velocity values, under the assumption that 




4820 



4840 



4860 
Wavelength (A) 



Figure 8. Comparison of the spectrum around H/3 observed at orbital phase 
0. 143 (blue points binned by two for clarity) to a composite synthetic spec- 
tram generated for the atmospheric parameters we found for the compo- 
nents of XY Cet and their expected radial velocities at this orbital phase 
(thin red line). 



the rotational axes of the stars are aligned with their orbital axis. 
The rotational velocity of star A is therefore synchronous (Tabled). 
StarB is rotating slightly faster than the synchronous velocity, al- 
though this is only significant at the 2.5<r level. 

The disentangled spectra of the stars contain their H7 lines, 
which are a useful indicator of T c g for A and F st ars but suffer fro m 
a degeneracy between T c g and surface gravity (Smalley 20 05bl iah. 
We have pursued two methods of obtaining T e g from the spectra, 
in both cases fixing the surface gravity (log g) values to those found 
from the masses and radii of the stars (see Sect. [7}, the relative light 
contributions of the stars to the value obtained from the light curve 
analysis (Sect.[5} and the v sin i values to the numbers above. 

Method 1 is a one-dimensional interpolation in T e g within 
a grid of LTE synthetic spectr a calculated using the SPECTRUIvQ 
code dGrav & Corballvl [l994h and ATLAS9 model atmospheres 
dKuruczll 19931) . A grid of spectra centred on H7 was constructed, 
covering T eff s from 7000 to 10 000 K in steps of 200 K. Residu- 
als were calculated between the observed and synthetic spectra, in 
carefully-selected spectral intervals which are free of metal lines. 
This approach yielded T c ft,A = 7870 ± 130 K and T*eff,B — 
7620 ± 170 K. 

Method 2 utilised the GENFITT code jTamaio et"al1 1201 ll) . 
which fits synthetic spectra to hydrogen line profiles using a ge- 
netic algorithm. This method yielded T c g A = 7870 ± 115 K and 
T e H,B = 7650 ± 125 K, in good agreement with method 1. A 
comparison between the observed and calculated H7 line profiles 
is shown in Fig.|7] An attempt to fit for log g as well as T e g re- 
turned less precise and somewhat different values, underlining a 
big advantage of dEBs that logy can be fixed in such analyses. 
We also obtained a solution with the light ratio fitted for and found 
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Figure 9. The observed and disentangled spectra of XYCet in the 4300-4340 A wavelength region. Star A is plotted using blue points offset by +0.4 and 
starB is plotted using green points. Synthetic spectra calculated for the appropriate T c ffS, log gs and V s i n iS and a solar chemical composition are shown by 
red lines. Some absorption lines are labelled with the species which causes them. 



0.72 ± 0.08, in agreement with the light curve results. The rela- 
tively large uncertainty arises because the H7 profiles contain many 
metal lines which have to be masked out in the GENFITT analysis, 
resulting in a significant loss of information. Better results could be 
obtained for the Ha or H/3 profiles, as these have far fewer metal 
lines, or for stars which are hotter or metal-poor. 

The T c ffS found above are in good agreement with those 
coming from the Stromgren indices (Table[7). We therefore adopt 
7870 ± 115 K and 7620 ± 125 K as our final values. The compo- 
nents of XY Cet turn out to be substantially cooler than was pre- 
viously thought. As a test we generated solar-abundance synthetic 
spectra covering the H/3 region, and combined them using the ap- 
propriate light ratio and RVs for our single observed spectrum of 
H/3. Fig.[8] shows that there is an excellent agreement between the 
synthetic composite spectrum and our observations. Our data rule 
out the hotter T ff s found in previous studies to high confidence. 



be found at 4283.011, 4318.562 and 4455.887 A, and Sell lines at 
4294.787, 4314.089, 4320.732, 4324.996 and 4400.389 A. These 
are uniformly weaker than solar, most obviously in the region of 
4320 A (Fig.[9j, so we confirm that both components are Am stars. 

Another intriguing facet of the Am phenom enon is a strong 
overabundance of rare earth elements (Wolff il983l) . We have there- 
fore carefully combed the disentangled spectra to uncover absorp- 
tion lines which are far stronger than expected for a solar com- 
position. A number of good examples have been found, compris- 
ing Lall 4333.753 A (see Fig.H, Cell (4364.653 and 4449.30 A) 
and Nd II (4446.384 A). Zr II (4403.349 A) is also highly enhanced. 
These elements have overabundances of order 1 dex in the photo- 
spheres of both components of XY Cet. This system is an excel- 
lent candidate for a detailed abundance analysis based on (echelle) 
spectra with a wide wavelength coverage and a high S/N ratio. 



6.1 Chemical abundances 

A major motivation for the current study was the metallic-lined 
nature of the components of XY Cet, which makes their physical 
properties useful in investigating the Am phenomenon. The limited 
wavelength coverage of our spectroscopic data unfortunately pre- 
cludes a detailed abundance analysis, but it is possible to obtain 
some qualitiative results by comparing the disentangled spectra to 
synthetic spectra calculated for a solar chemical composition. We 
adopted a microturburbulence velocity of 2 km s ~ 1 but our conclu- 
sions are robust against a n increase to a value of 4km s~ 1 , as often 
found for Am stars (e.g. ICouprv & Burkhadll992l Fig.|9] shows 
this comparison for the 4300-4340 A wavelength region, with a se- 
lection of absorption lines labelled. 

From Fig.|9]it is obvious that most of the spectral lines are 
significantly enhanced compared to solar abundance, which is a 
common feature of Am stars. The prime hallmark of the Am phe- 
nomenon is that Ca and Sc are underabundant in the photosphere, 
causing their absoiption lines to be unusually weak. Ca I lines can 



7 THE PHYSICAL PROPERTIES OF XY CET 

We have calculated the physical properties of the XYCet sys- 
tem based on the results from the photometric and spectro- 
scopic analyses above. T his was done with the JKTABSDIM code 
( Southwort h et al.l2005al) , which calculates complete error budgets 
using a perturbation algorithm. The physical properties of XY Cet 
are given in Tableland show that we have measured the masses 
of the stars to accuracies of 1% and the radii to accuracies of 2%. 
Additional photometric observations could allow the radius mea- 
surements to be improved to the 1% level or better. A spectroscopic 
light ratio would be beneficial, and could most usefully be obtained 
from disentangling spectra covering the H/3 line. 

Using the T ff s measured in Sect.|6]we have measured the dis- 
tance to the XY Cet system using the method of ISouthworth et all 
(2005a), which is based on calibrations be tween T e ff and surface 
brightness obtained from interferometry dKervella et al.ll2004h . An 
interstellar reddening of Eb-v = 0.04 ± 0.02 is needed to obtain 
agreement between the distances found in the optical (UBV) and 
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Table 8. The physical properties of the XY Cet system. 



Parameter 


Star A 


StarB 


Orbital separation ( Rq) 


12.493 ± 0.035 


Mass (Mq) 


1.773 ±0.016 


1.615 ±0.014 


Radius ( Rq ) 


1.873 ±0.035 


1.773 ±0.029 


iu & y L Lm h j 


4.142±0.016 


4.149 ±0.014 


V synch (kms" 1 ) 


34.09 ±0.64 


32.27 ±0.53 


V s im (kms -1 ) 


34.4 ±0.4 


34.1 ±0.5 


T cff (K) 


7870 ±115 


7620 ±125 


log(L/L Q ) * 


1.082 ±0.030 


0.978 ±0.032 


M hoi * 


2.045 ± 0.076 


2.304 ±0.080 


Distance (V-band, pc) 


285 ±13 


292 ± 14 


Distance (A"-band, pc) 


279.3 ± 5.6 



* Calculated ass uming hp, = 3.844xl0 26 W jBahcall et aj]|l995h and 
M bo i0 = 4.75 <Zombeckll99"rJI) . 

near-infrared (JHK) passbands. The most precise distance mea- 
surement is found in the if-band, primarily because of its relative 
insensitivity to reddening and T e «, and is 279 ± 6pc. Consistent 
distances are found using the bolometric-correction technique (e.g. 
ISouthworth et al.l 2005b) and empirical or theoretica l bolo metric 
corrections from ICode et al.l dl976f) and iBessell et al.l (119981) . The 
good agreement between the distances for the two stars in the same 
passbands is supporting evidence that the light ratio obtained from 
the light curve is correct. 

A comparison of the masses, radii and T c ff s of X Y Cet to sev- 
eral sets of theoretical models of stellar evolution dClarellll995l 
l2004bl : iGirardi et alj|2000t IPols et al.|[l997h has been performed. 
We find that the best match corresponds to a solar helium and metal 
abundance and an age in the region of 850 ± 100 Myr, but the mass- 
radius relation of the dEB is less steep than predicted by any of the 
models. An enhanced helium or metal abundance can improve the 
agreement in the mass-radius plane, but at the expense of the agree- 
ment in the mass-T c fr plane. We conclude that the components of 
XY Cet appear to be discrepant with respect to theoretical predic- 
tions, but that the discord is smallest for solar-composition models. 



8 SUMMARY 

XYCet is an eclipsing binary containing two components which 
show the Am phenomenon. It has been studied many times previ- 
ously, but without agreement on the radii of the stars and based 
predominantly on photographic spectra. We have measured the 
physical properties of the two stars from new phase-resolved CCD 
spectroscopy and an extensive light curve obtained by the Super- 
WASP survey for transiting planets. The spectra were analysed us- 
ing cross-correlation and spectral disentangling techniques, and the 
light curve was modelled using the JKTEBOP code. We measured 
the effective temperatures of the stars by fitting their H7 line pro- 
files, finding good agreement with values from calibration against 
Stromgren indices. 

Our spectroscopic data cover only a narrow wavelength re- 
gion so are not useful for a detailed chemical abundance analysis. 
Qualitative conclusions have therefore been obtained by comparing 
the disentangled spectra of the two components to synthetic spec- 
tra calculated for the atmospheric parameters of the stars. We find 
that both stars demonstrate a clear overabundance of iron-group 
elements and an obvious deficiency of calcium and scandium, as 
expected for Am stars. Further confirmation of their metallic-lined 



nature comes from the detection of strongly enhanced abundances 
of zirconium and the rare-earth elements lanthanium, cerium and 
neodymium. 

The physical properties of XY Cet are measured to accuracies 
of 1% in mass, 2% in radius and 125 K in T c a ■. We find that the stars 
are rotating close to synchronicity with the orbital velocity. The 
system is located at a distance of 279 ± 6 pc. Theoretical models 
of stellar evolution are able to match the T c fts of the stars for the 
measured masses, but predict a steeper slope in the mass-radius 
diagram than observed. 

The physical properties of the component stars of XY Cet are 
not clearly different from those of normal A stars, supporting the in- 
terpretation of the Am phenomenon as a surface effect which does 
not significantly affect the evolution of stars. The chemical pecu- 
liarities and relatively slow rotation of both components mean that 
XY Cet is a promising target for an extensive spectroscopic inves- 
tigation of the photospheres of metallic-lined stars with accurately 
known masses, radii and T c s s. 
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